Stability of inhomogeneous superstructures from renormalized mean-field theory of 
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Using the t— J model (which can also include Coulomb repulsion) and the "plain vanilla" renor- 
malized mean-field theory of Zhang et al. (1988), stability of inhomogeneous 4a x 4a superstructures 
as those observed in cuprates superconductors around hole doping 1/8 is investigated. We find a 
non-uniform 4a x 4a bond order wave involving simultaneously small (~ 10 -2 t) inhomogeneous 
staggered plaquette currents as well as a small charge density modulation similar to pair density 
wave order. On the other hand, no supersolid phase involving a decoupling in the superconducting 
particle-particle channel is found. 

PACS numbers: 75.10.-b, 75.10.Jm, 75.40.Mg 
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Although the global phase diagram of high-Tc cuprate 
superconductors seems, at first sight, to contain only 
three homogeneous phases, antiferromagnetic, supercon- 
ducting and metallic phases, a closer inspection reveal a 
striking complexity with some form of local electronic or- 
dering around 1/8 hole doping. For exemple, recent scan- 
ning tunneling microscopy/spectroscopy (STM/STS) ex- 
periments of under-doped Bi2Sr2CaCu208+<5 in the pseu- 
dogap state have shown evidence of real-space modula- 
tions of the low-energy density of states (DOS) Q with 
period close to four lattice spacings. Spatial variation 
of the electronic states has also been observed in the 
pseudogap phase of C^-zNazCuC^C^ single crystals 
(x = 0.08 - 0.12) by similar STM/STS techniques 0. 
In the same material, signatures of this unusual 4a x 4a 
checkerboard charge-ordered state have been found in 
Angular Resolved Photocmission Spectroscopy (ARPES) 
experiments How this local ordering could be ex- 
plained on theoretical grounds is a topic of intense ongo- 
ing investigation. 

In the Resonating Valence Bond (RVB) theory Q, 
one of the early approach to understand strongly corre- 
lated HTC superconductors, a simple variational Ansatz 
is constructed from a BCS paired superconducting state 
I^bcs) by a Gutzwiller [5J projection of the high energy 
configurations with doubly occupied sites. At half-filling, 
this procedure generates a highly correlated wavefunc- 
tion PgI^bcs) where the spins are paired up in singlet 
bonds to form a quantum spin liquid with short range 
antiferromagnetic correlations. In fact, due to a SU(2) 
electron-hole symmetry, the RVB state can take various 
forms like the half-flux state which can be mapped 
onto free electrons on a lattice experiencing half a flux 
quantum per plaquette. 

The main difficulty to deal with such wavefunctions 
is to treat correctly the Gutzwiller projection. Exact 
Variational Monte Carlo (VMC) calculations [jj on large 
clusters show that the magnetic energy of the variational 
RVB state at half-filling is very close to the best ex- 



act estimate. It also provides, at finite doping, a semi- 
quantitative understanding of the phase diagram of the 
cuprate superconductors and of their experimental prop- 
erties. Another route to deal with the Gutwillcr pro- 
jection is to use a "renormalized mean-field (MF) the- 
ory" in which the kinetic and superexchange energies 
are renormalized by different doping-dependent factors gt 
and gj respectively. As emphasized recently by Anderson 
and coworkers Q , the procedure consists in searching for 
an effective MF Hamiltonian that determines a function 
I*) to be projected as Pg\^). Crucial, now well estab- 
lished, experimental observations such as the existence 
of a pseudo-gap and nodal quasi-particles and the large 
renormalization of the Drude weight are remarkably well 
explained by this early MF theory [jj. 

Away from half-filling, "commensurate flux phases" 
have also been proposed as another class of competi- 
tive projected wavefunctions 0,0]. Such time rever- 
sal symmetry broken states, generalizing Affleck-Marston 
1/2-fhix phase away from half- filling, are made out of 
orbitals solving a fictitious Hofstadter problem 0] of 
free particles under a (fictitious) magnetic flux. Within 
this class of wavefunctions, the lowest energy is obtained 
when the flux per unit cell equals exactly the filling frac- 
tion v = |(1 — x) (x being the hole doping) 0, 
suggesting that doping induces a (scalar) chirality of the 
spin background [1J|. This scenario also emerges nat- 
urally within a renormalized MF theory |l5j although 
the stability condition (self-consistency) leads to a more 
complex one-body effective MF Hamiltonian with inho- 
mogeneous flux. For example, at filling v = p/q, diagonal 
"stripe" modulations of supercell qa/^/2 x \J~2a [a is the 
lattice spacing) are found. In this Letter, we generalize 
the MF approach of Ref. to allow for non-uniform 
densities (and also include off-diagonal pairing). We de- 
termine under which conditions a 4a x 4a superstructure 
might be stable for hole doping close to 1/8. 

The weakly doped antiferromagnet is described here 
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by the renormalized t— J model Hamiltonian, 

H = -tgt J2 (4,*Cj,« + + J9J E S < ' S - (!) 

where the local constraints of no doubly occupied sites 
are replaced by statistical Gutzwiller weights gt — 
2x/(l + x) and gj — 4/(1 + a;) 2 . We shall assume a 
typical value of t/J = 3 hereafter. 

Decoupling in both particle-hole and (singlet) particle- 
particle channels can be considered simultaneously lead- 
ing to the following MF hamiltonian, 

Hmf = ~* E 4K c i,ff c J> + /l - c ') ~~ f x ^2 n h" 

icr 

~\ J E 9ij(Xji4 ta c jt(7 + h.c. - \ Xt] \ 2 ) (2) 

-\ J E gl^AA,-* + h - c - - \ A v I 2 )' 

where the previous Gutzwiller weights have been ex- 
pressed in terms of local fugacities Zi — 2xi/(l + SEj) (xi 
is the local hole density), gjj = ^Jziz] and gf, = (2 — 
Zj)(2 — Z i), to allow for small non-uniform charge mod- 
ulations [lfij . The Bogolubov-de Gennes self-consistency 
conditions are implemented as Xji — { c \,<j Ci > a ) anc ^ = 

Assuming an homogeneous system, the RVB theory of 
Refs. 0,01 is recovered leading to a d x i_ y i RVB super- 
conducting (for x > 0) ground state (GS), A; i +x = A, 
A, i+ y = —A and Xij = X- When doping x — * 0, 
A = x — 0.169, giving a magnetic energy per bond 
of ~ 0.339J and a kinetic energy per hole of ~ -2.7U 
(within MF). J A is physically the pseudo-gap while the 
superconducting temperature scales like gt J A oc x. 

Our aim here is to investigate inhomogeneous solutions 
of these equations which might a priori break translation 
symmetry. However, searching for the global energy min- 
imum in a completely unconstrained variational space is 
a formidable task owing to the non-linear nature of the 
MF equations. Therefore, guided by experimental ob- 
servations, for an average hole density of x ~ 1/8, we 
shall restrict ourselves to a 4a x 4a superstructure (which 
could, as well, accommodate a smaller spatial periodic- 
ity) with a natural C4 rotation symmetry, hence natu- 
rally restricting the number of non-equivalent bonds and 
sites to 6 and 3, respectively. Within these restrictions, 
MF equations are solved on finite clusters (with periodic 
boundaries) of sizes ranging from 16 x 16 to 48 x 48. 

We first consider non-superconducting solutions i.e. 
Aij = on all bonds. Starting with an initial set of 
Xij which simultaneously (i) corresponds to a uniform 
(fictitious) flux (per plaquette) equal to the filling frac- 
tion v = 7/16 and (ii) fulfills the above symmetry re- 
quirements within the supercell, we then solve the self- 





FIG. 1: (Color on-line). Schematic pattern of the bond Xij 
distribution. Non-equivalent bonds (labelled from 1 to 6) are 
indicated by different types of lines whose widths qualitatively 
reflect the magnitudes of Arrows indicate the directions 
of the small charge currents. Numerical values for t/J = 3 
are given in Table [I] 

consistent equations iteratively 01- A typical MF solu- 
tion is shown schematically in Fig. ^ As shown in Ta- 
ble |U the bond 1 mramctcrs Xiji here obtained e.g. on a 
48 x 48 cluster |l8j , show strong spatial modulations lead- 
ing to a bond order wave (BOW) in both spin-spin cor- 
relations fpylxyl 2 an d charge hoppings 2g t i] Real{xij}- 
At first sight, such a state could be seen as an extension 
of the bond-centered stripe 0] or 4 x 4 plaquette 
solutions found in some SU (2N) /Sp(2N) mean field the- 
ories. However, in our case, xt.j exhibit imaginary parts 
(as in staggered |2l| and commensurate flux states) pro- 
ducing a time-reversal symmetry broken state and small 
plaquette charge currents 2i?* J Im{xr,} (see Table|IJ. The 
current pattern is similar to the alternating distribution 
of the staggered flux state although the amplitude of the 
current loops is not uniform. As shown on the first line of 
Table UTI such a state also exhibits a small charge density 
wave (CDW) component. 

TABLE I: Numerical values of the bond amplitudes |x*i I> the 
bond spin-spin correlations and the bond average hoppings 
and currents (in units of t) in the normal phase MF solution 
obtained on a 48 x 48 cluster for t/J = 3 and 2022 fermions 
on 48 2 sites (x ~ 1/8). Note the current conservation at the 
lattice sites, in particular "2+6" ="3+4" . 



Bond # 


1 


2 


3 


4 


5 


6 


\xn\ 


0.2579 


0.1871 


0.2595 


0.2259 


0.1887 


0.1639 


(Si -Sj) 


-0.3456 


-0.1720 


-0.3131 


-0.2359 


-0.1637 


-0.1249 


hoppings 


0.0554 


0.0571 


0.1188 


0.1063 


0.0915 


0.0725 


currents 


0.0453 


0.0373 


0.0370 


0.0305 


0.0228 


0.0302 
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FIG. 2: (Color on-line). Schematic pattern of the charge and 
bond modulation. Non-equivalent sites (labeled from 1 to 3) 
and bonds are indicated by different types of circles and lines 
(and colors, online). Numerical values for t/J = 3 are given 
in Table El 

We now turn to superconducting solutions allowing 
pair order parameter Ay ^ 0. We repeat the previ- 
ous procedure adding to Xij arbitrary initial values for 
Ay, only requesting |Ay| to have s-wave symmetry w.r.t 
the center of the central plaquette (defined by bonds #1). 
Within this restriction, we have looked for 3 candidates 
compatible with a 4a x 4a cell, exhibiting (i) d x 2_ y 2 +is, 
(ii) d x 2_ y -2 + id xy or (iii) d x 2_ y 2 + id x 2_ y 2 symmetries. 
Note that we allow again inhomogeneous solutions so 
that the minimization of the MF energy w.r.t. the set 
of Ay is performed indeed on 6 independent complex 
parameters. However, for a doping around 1/8 and all 
cluster studied, no modulated solution was found. On 
the contrary, the iteration procedure always converges to- 
wards the standard d x 2_ y 2 RVB superconductor, a state 
particularly rigid w.r.t. the establishment of spontaneous 
modulations. 

Next, we investigate the role of a Coulomb repulsion, 
H v = ^ Vl >J ( ni ~ n ) ( n i ~ n ) ' ( 3 ) 

that we treat within a simple decoupling scheme, i.e. 
replacing the density-density interaction by (ni)rij + 
(nj)ni + Cst. We assume here both long range Vij = 
V/rij and screened V%j = V exp (— ry /0/ r »,j potentials 
(with a screening length set to I — 4a). To minimize fi- 
nite size effects, we also use the "periodic distance" ry 
on the torus. As shown in Table ITT1 we find that moder- 
ate values of V (or V') ~ t affects only slightly the form 
of the CDW. Moreover, the pattern for both spin-spin 
correlations and effective hoppings remains qualitatively 
similar to V — (see Table HID . Note that the charge 
current magnitude tends to increase (decrease) with V 



on the central, (outer) plaquette, i.e. on bonds #1 (#5). 
As expected, results do not depend drastically on the 
long-distance part of Vij . 

TABLE II: Numerical values of charge densities in the normal 
phase obtained on a 48 x 48 cluster for t/J = 3 and doping 
x ~ 1/8 and different Coulomb repulsions. Average density, 
magnetic bond energy, kinetic energy per hole and energy per 
site are provided in the 4 last columns. Last line: same data 
obtained for the d-wave RVB superconductor on a 32 x 32 
cluster for fj, = -0.4 (leading to X ^ 0.192 and A ~ 0.120). 



Site # 


1 


2 


3 


n 


(Si ■ 8j) 


E K /t 


e /J 


V = 


0.925 


0.864 


0.857 


0.8776 


-0.220 


-2.660 


-1.438 


V/t = 0.5 


0.894 


0.913 


0.783 


0.8759 


-0.227 


-2.588 


-1.425 


V/t = 0.8 


0.876 


0.940 


0.760 


0.8793 


-0.232 


-2.396 


-1.367 


V'/t = 0.8 


0.850 


0.932 


0.786 


0.875 


-0.222 


-2.559 


-1.406 


RVB 


0.876 


0.876 


0.876 


0.8762 


-0.243 


-2.730 


-1.500 



Let us now turn to the issue of the relative stability 
of the BOW/CDW state w.r.t. the RVB superconduc- 
tor. For V — the latter has a slightly lower energy. 
However, a moderate repulsion can easily stabilize the 
BOW/CDW state whose energy is only weakly affected 
by V (or V). On the contrary, the pairing amplitude of 
the RVB superconductor is expected to be quite sensitive 
to the Coulomb potential. Indeed, a decoupling in the 
particle-particle channel of e.g. the nearest neighbor re- 
pulsion Vin^Uji term generates a positive contribution 
tx Vi|Ay| which tends to strongly reduce the pairing 
amplitude. Therefore, the BOW/CDW state is expected 
to become stable when x ~ 1/8. 



TABLE III: Same as Table U] but for a screened Coulomb 
potential V'/t = 0.8 (/ = 4a). 2016 fermions on 48 2 sites 
(x = 1/8). 



Bond # 


1 


2 


3 


4 


5 


6 


IXijI 


0.2297 


0.2015 


0.2741 


0.1917 


0.2166 


0.2020 




-0.2395 


-0.1984 


-0.3952 


-0.1700 


-0.1909 


-0.2146 


hoppings 


0.1105 


0.0637 


0.0547 


0.0773 


0.1524 


0.0403 


currents 


0.0459 


0.0366 


0.0436 


0.0251 


0.0109 


0.0321 



Lastly, we would like to discuss similarities and differ- 
ences of our findings with previous theoretical suggestions 
to explain the 4a x 4a patterns experimentally observed 
e.g. in BiSCO. First, we should stress that the state 
found in this work is quite different from an ordinary 
CDW or Wigner solid which, for a doping around 1/8, 
is expected to have a x v° supercell. As a matter 
of fact, the BOW is not primarily stabilized by Coulomb 
repulsion. Secondly, we notice that it bears some simi- 
larities (in addition to the common 4a x 4a periodicity) 
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with the insulating Pair density Wave (PDW) state [2^ 
which also does not exhibit U(l) phase coherence. For 
exemple, it shows a higher hole density on a given pla- 
quette within the supcrcell (which, incidcntly, carries the 
minimum value of the current). Like the PWD state, the 
BOW differs from the stripe state in terms of the rota- 
tional symmetry. In addition, in contrast to PDW and 
stripe states, it breaks time-reversal symmetry. However, 
an experimental observation of its manifestation i.e. the 
plaquette orbital currents might be quite tedious due to 
their small magnitudes. We note that the current pat- 
tern is, at first glance, staggered (in fact, incommensurate 
with a wave vector close to (7r,7r)), a finding consistent 
with exact diagonalisations of doped Hubbard and t— J 
clusters showing enhancement of staggered chiral correla- 
tions under doping |23|. Finally, we observe that recent 
VMC calculations of staggered flux states [U suggest 
an instability towards phase separation. A modulated 
structure like the one found here is in fact an alternative 
scenario. 

In conclusion, using renormalized mcanfield theory of 
the t— J model, 4a x 4a BOW states are found for dop- 
ing around 1/8. We suggest that such states might be 
relevant to explain the spatial modulations of the DOS 
in STM experiments of underdoped cuprates supercon- 
ductors. They carry small (almost staggered) plaquette 
currents and are stable for moderate Coulomb repulsion. 
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